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INTRODUCTION 


GLOBALL  was  developed  to  give  the  desk-bound  navigator  computer-aided  solutions 
to  his  global  range  and  surface  area  problems.  The  program  can  accomplish  tasks  that 
are  tedious  and  time  consuming,  if  not  virtually  impossible,  to  perform  by  hand.  To 
verify  inputs  and  visualize  results,  a  three-dimensional,  computer -drawn  plot  of  the 
earth’s  surface  is  furnished.  The  plot  can  also  be  used  for  graphic  presentation  of  work 
in  a  document. 

Three  basic  calculations  can  be  performed:  determining  great-circle  distances,  gen¬ 
erating  range  circles  about  points,  and  computing  surface  areas  from  inputted  point  coor¬ 
dinates  describing  the  perimeters  of  the  arqas.  These  operations  can  be  combined  to 
solve  a  variety  of  problems . 

GLOBALL  was  designed  to  be  used  as  a  library  routine  to  call  various  subroutines 
which  use  inputs  that  are  read  in  and  stored,  or  values  that  are  calculated  during  a  run. 

This  paper  presents  a  general  description  of  the  basic  computer  program  and  its 
accompanying  plotter  subroutine,  a  detailed  description  of  the  inputs  required  in  the  sub¬ 
routines,  and  several  examples  of  possible  applications.  Each  example  is  accompanied  by 
a  program  listing,  computer  printouts,  and  illustrative  plots. 


PART  I:  GENERAL  DESCRIPTION 

The  GLOBA LL  Computer  Program,  which  is  a  package  of  Fortran  IV  subroutines,  was 
developed  at  the  Center  for  Naval  Analyses  to  compute  global  distances  and  areas  anywhere 
on  the  surface  of  the  earth.  It  evolved  from  three  computerized  algorithms  that  were  used 
for  calculating  great -circle  distances,  range  circles,  and  surface  areas.  The  algorithms 
were  combined  with  a  specially  adapted  three-dimensional  plotting  subroutine,  forming  the 
three  main  subroutines:  GCLINE,  SURFAREA,  and  GLOBEMAP. 

Great-Circle  Distances 

SUBROUTINE  GCLINE  calculates  great -circle  distances  between  points  on  the  earth's 
surface  in  one  degree  segments  (60  n.mi.)  and  plots  them.  The  coordinates  at  the  junc¬ 
tion  of  each  segment  are  computed  and  printed  out. 

Range  Circles 

A  range  circle  may  be  generated  about  any  point  on  the  earth's  surface  and  is  formed 
by  revolving  its  great -circle  radius  about  its  center.  The  track  of  the  range  circle  is 
divided  into  one  degree  arcs  beginning  at  true  north  clockwise  through  360  degrees.  A 
number  of  range  circles  may  be  calculated  and  drawn  on  the  same  plot.  If  the  circles 
intersect,  the  intercept  points  and  the  areas  of  intersection  are  computed. 

Surface  Areas 

SUBROUTINE  SURFAREA  will  compute  a  surface  area  of  the  portion  of  the  earth's 
surface  enclosed  within  a  given  perimeter.  A  perimeter  described  by  a  set  of  points  may 
represent  an  area  up  to  a  hemisphere  in  size.  An  area  enclosure  may  take  on  any  shape 
and  cover  both  land  and  water. 

Projection  Map 

GLOBALL  provides  an  accurate  and  easy-to-use  analytical  tool  for  use  in  a  wide 
variety  of  applications.  To  get  a  better  "feel"  for  the  numerous  possible  range  and  area 
configurations,  a  two-dimensional,  global  representation  or  map  is  plotted  by  the 
SUBROUTINE  GLOBEMAP. 

GLOBEMAP  PLOTTER  SUBROUTINE 

GLOBEMAP  was  developed  from  another  three-dimensional  plotting  subroutine  that 
was  used  to  plot  quadric  surfaces,  and  draws  azimuthal  orthographic  projections  of  the 
earth's  curved  surface.  Either  a  full  hemisphere  or  an  enlarged  projection  may  be  selected. 
The  physical  size  of  the  plots  may  be  varied.  The  projection  plots,  showing  the  coastlines 
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of  the  earth  land  masses,  make  excellent  document  maps  with  titles  along  three  of  the 
four  sides. 

Inheriting  the  characteristics  of  an  azimuthal  projection,  every  point  is  shown  at  its 
true  direction  from  the  center  of  the  projection.  Similarly,  a  great -circle  distance  from 
the  center  point  to  any  other  point  projects  as  a  straight  line. 

The  user  has  complete  freedom  to  specify  any  point  on  the  surface  of  the  earth  as  the 
center  of  the  projection.  Each  view  is  obtained  after  rotating  the  globe  about  fixed  spatial 
axes. 

Spherical  coordinates  (p,  6  ,  <t>)  are  used  to  pinpoint  the  positions  of  points  lying  on  the 
projected  sphere.  The  radius  of  the  sphere,  p,  is  assumed  to  be  one  unit  in  length  to 
simplify  calculations.  To  compute  surface  areas,  P  is  set  equal  to  the  earth’s  radius, 
3440.07  nautical  miles.  0and  4>  are  the  co-latitude  and  longitude  of  a  point,  respectively. 

Types  of  Plots 


The  user  has  the  option  to  draw  either  a  globe  projection,  showing  a  hemisphere  at  a 
time,  or  an  enlarged  view  of  a  smaller  part  of  the  earth’s  surface.  Examples  of  both  are 
illustrated  in  part  III  of  this  paper. 

The  hemisphere  projection  is  drawn  by  plotting  only  the  visible  half  of  the  earth. 

Plotting  only  the  points  that  are  within  5,  400  nautical  miles  of  the  center  of  the  plot,  the 
back  side  is  suppressed.  A  circular  border  or  horizon  is  drawn  with  a  radius  of  5,  400  n.mi. 
SUBROUTINE  GLOBEMAP  is  called  for  the  hemisphere  plot. 

A  separate  subroutine,  BLOWUP,  is  called  for  enlarged  projections.  An  input  param¬ 
eter,  BRANGE,  which  will  be  referred  to  as  the  blow-up  range,  enables  the  user  to  vary 
the  magnification  of  the  projection.  The  blow-up  range  is  the  maximum  distance  from  the 
center  of  the  projection  outside  of  which  points  are  not  plotted.  This  distance  is  internally 
set  equal  to  5,  400  n.mi.  for  a  hemisphere  plot.  The  reduced  field  of  view  is  then  enlarged 
to  completely  fill  a  circular  or  rectangular  border.  The  dimensions  of  the  border  are 
preset  by  other  inputs  and  unchangeable  during  the  enlarging. 

The  ability  to  vary  the  magnification  and  the  physical  size  of  a  projection  enables  the 
user  to  "zoom"  in  on  a  particular  region. 

Shorelines 


The  earth's  land  masses  and  inland  bodies  of  water  are  drawn  by  connecting  shoreline 
points.  The  mean  distance  between  the  points  is  approximately  8  nautical  miles.  All  islands, 
lakes  and  inland  seas  greater  thanfc  250  square  nautical  miles  are  represented. 
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The  coordinates  of  the  points  are  stored  on  tape  and  read  in  for  each  plot.  Each  point 
of  every  shoreline  is  tested  to  determine  whether  it  is  within  the  dimensions  of  the  plot.  A 
shoreline  is  drawn  from  point  to  point  until  either  all  its  points  have  been  read  in  or  until 
it  extends  past  the  border.  A  shoreline  may  enter  and  exit  a  plot  several  times. 

Meridians  and  parallels,  spaced  every  ten  or  twenty  degrees,  may  be  added. 
CONCLUDING  REMARKS 

The  ability  to  rotate  the  projections  of  the  earth's  surface  and  draw  a  sequence  of 
plots  adds  a  fourth  dimension  --  time.  Production  of  time -variable  plots  enables  the 
researcher  to  pictorially  analyze  campaign  results  at  various  stages. 

A  future  version  of  GLOBALL  is  currently  under  development  to  improve  the  con¬ 
ceptual  quality  of  the  projection  maps  by  including  the  political  boundaries  of  the  countries. 
Two  additional  subroutines,  RADARC  and  FIRANGE,  will  be  added  to  compute  and  draw 
radar  azimuthal  coverages  and  missile  firing  ranges,  respectively. 

GLOBALL  can  provide  a  valuable  graphic  tool  to  visualize  and  study  satellite  ground 
track  patterns.  The  satellite's  path  across  the  surface  of  the  earth  can  be  plotted  by 
SUBROUTINE  GCLINE  and  coverage  areas  computed  by  SUBROUTINE  SURFAREA.  Addi¬ 
tional  modifications  to  GLOBA  LL's  coordinate  transformation  routine  would  be  required  to 
actually  plot  the.  satellite's  orbit  and  any  other  space  trajectories. 
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PART  II:  INPUT  PARAMETERS  AND  SUBROUTINES 


INPUTS 

The  principal  inputs  are  coordinates  of  points,  (<r,  \  )  and  ranges  (n.mi.).  The  coor¬ 
dinates  identify  the  positions  of  points  for  which  an  indicated  operation  is  to  be  performed. 
The  latitudinal  position  of  a  point  measured  in  degrees  is  the  angular  distance  north  or 
south  of  the  equator.  As  an  input  parameter,  the  latitude,  <r,  of  a  point  is  considered 
positive  north  of  the  equator  and  negative,  south.  Similarly,  longitude,  X,  is  measured 
positive,  east,  of  the  prime  meridian. 

GLOMLL's  subroutines  and  input  parameters  are  listed  as  follows: 

GLOBEMAP  (XLAT,  XLON,  PRANGE,  NOPTS,  XTITLE,  NX,  YTITLE,  NY, 

TITLE,  NT,  NTYPE,  XSIZE,  YSIZE) 

BLOWUP  (XLAT,  XLON,  PRANGE,  BRANGE,  NOPTS,  XTITLE,  NX,  YTITLE, 

NY,  TTITLE,  NT,  NTYPE,  XSIZE,  YSIZE) 

SECRANGE  (XLAT,  XLON,  SRANGE,  NOPTS) 

GCLINE  (XLATS,  XLONS,  NPTS,  NOPTS) 

SURFAREA  (XLATS,  XLONS,  NPTS,  CX,  CY,  NOPTS)1 

The  list  of  inputs  consists  mostly  of  the  usual  plotting  parameters.  XTITLE,  YTITLE, 
and  TTITLE  are  used  for  placing  titles  containing  NX,  NY  and  NT  number  of  characters 
along  three  of  the  four  sides  of  a  plot.  NTYPE  orients  the  plot  for  viewing  from  either 
side  of  the  plot  page.  For  NTYPE  =  0,  the  plot  is  viewed  from  the  11"  dimension;  for 
NTYPE  =  1  from  the  8-^"  dimension.  XSIZE  and  YSIZE  are  the  plot's  dimensions  expressed 
in  inches. 

NOPTS  is  an  array  of  inhouse  options  built  into  GLOBA LL  to  display  internal  computa¬ 
tions,  alter  or  suppress  printouts,  and  control  various  aspects  of  the  projection  maps. 

None  of  the  options  were  employed  for  the  examples  in  this  paper. 


SURFAREA's  inputs  are  listed  and  described  on  page  31. 
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Other  inputs  are: 


•  GLOBEMAP 

XLAT,  XLON  =  latitude  and  longitude  of  center  point  of  plot  and  center 
of  a  primary  range  circle. 

PRANGE  =  radius  of  a  primary  range  circle  in  nautical  miles 

•  BLOWUP 

XLAT,  XLON  =  same  as  above. 

PRANGE  =  same  as  above. 

BRANGE  =  blow-up  range  (discussed  later  in  the  paper). 

•  SECRANGE 

XLAT,  XLON  =  center  coordinates  of  secondary  range  circle. 

SRANGE  =  radius  of  secondary  range  circle  in  nautical  miles. 

•  GCLINE 

XLATS  =  array  of  latitudes  of  NPTS  number  of  points 
XLONS  =  array  of  longitudes  of  NPTS  number  of  points 
NPTS  =  number  of  points 

AUXILIARY  FUNCTION  SUBROUTINES 

All  spherical  distances  and  angles  used  by  the  previous  mentioned  subroutines  are 
computed  by  the  function  subprograms  CALCSIDE,  CALCANGL  and  OPPOSIDE.  FUNCTIONS 
CALCSIDE  and  CALCANGL  are  used  together  to  determine  great -circle  distances  and  range 
circles.  FUNCTION  OPPOSIDE  is  the  spherical  trigonometric  identity  required  for  com¬ 
puting  intercept  points  of  intersecting  ranges. 

Given  two  sides  and  an  included  angle  of  a  spherical  triangle,  FUNCTION  CALCSIDE 
calculates  the  third  side . 

CALCSIDE  =  COS"1  [(COS(Xl)  *  C0S(X2)  +  SIN (XI)  *  SIN(X2)  •  COS(X3)]  , 
where  XI  and  X2  are  the  two  sides  and  X3  is  the  included  angle. 


-6- 


Given  three  sides  of  a  spherical  triangle,  FUNCTION  CALCANGL  calculates  an  angle 
opposite  one  of  the  sides. 


CALCANGL  =  COS 


COS(X3)  -  COS(Xl)  ‘  COS(X2) 
SIN (XI)  *  SIN(X2) 


where  XI,  X2,  and  X3  are  the  three  sides  and  the  computed  angle  is  opposite  side  X3  . 


Given  two  sides  and  an  opposite  angle  of  a  spherical  triangle,  FUNCTION  OPPOSIDE 
calculates  the  third  side. 


OPPOSIDE  =  SIN’1  [SIN(TAN_1(TAN(X3)  •  COS  (XI)))  •  TAN(X1)  •  COT(X2)] 
-TAN"1  [TAN(X3)  •  COS(Xl)]  , 
where  XI  and  X2  equal  the  two  sides  and  X3  is  the  opposite  angle. 

MAP  CHARACTERISTICS 

To  draw  a  two-dimensional  representation  of  the  earth's  curved  surface,  all  points 
defined  in  spherical  coordinates  (P,  0,  <t>)  must  be  transformed  to  two-dimensional  carte¬ 
sian  coordinates  (X,  Y).  The  two  coordinates  systems  are  related  through  the  equations 

X  =  p-  SIN(0)  •  SIN(<*>) 

Y  =  P-  [  COS  (a)  •  COS(0)  -  SIN(a)  •  COS(<f>)  •  SIN(0)]  , 
where  a  =  latitude  of  the  center  point  of  the  projection. 

The  power  of  magnification,  X^,  of  an  enlarged  projection  is  inversely  proportional  to 
the  sine  of  the  degree  equivalent  of  the  blow-up  range. 


1.0 

SIN  (/3) 


where  0  =  BLOW-UP  range/60.0  (n.mi.  per  degree). 


Whether  a  circular  border  or  a  rectangular  border  is  drawn  around  a  map  is  deter¬ 
mined  by  the  blow-up  range.  For  the  rectangular  border,  this  range  is  measured  along 
the  longer  dimension.  The  range  must  be  small  enough  to  insure  that  the  diagonal  dis¬ 
tance  is  less  than  or  equal  to  4,  800  nautical  miles.  The  horizon,  which  is  5,  400  n.mi. 
from  the  center  of  the  plot,  is  not  allowed  to  appear  in  the  comers  of  the  plot.  A  600  n.mi. 


margin  is  observed  for  interpolation  to  determine  intercept  points  with  the  border.  If  the 
user  is  in  doubt,  the  following  hand  calculation  should  be  made  before  submitting  a  run. 


DlA=^z , 

where 

DL  =  longer  dimension, 

DS  =  shorter  dimension. 

Then, 

BLOW-UP  range  s  SIN  *(DIA  *  SIN(80°))  •  60.0  (n.mi.  per  degree). 

Example:  For  an  8"  x  6"  border,  the  blow-up  range  must  be  less  than  or  equal  to  3, 119 
nautical  miles. 


-8- 


PART  III:  APPLICATIONS  AND  EXAMPLES 


GREAT  CIRCLE  DISTANCES 

Great-circle  distances  can  be  computed  by  calling  the  SUBROUTINE  GCLINE.  The 
original  version  of  GCLINE  was  used  to  compute  a  great-circle  distance  in  nautical  miles 
between  two  points  on  the  earth's  surface.  Before  inclusion  as  a  subroutine  of  GLOBALL, 
GCLINE  was  given  the  additional  capability  to  be  called  with  two  or  more  points.  For  N 
number  of  points,  N-l  great-circle  distances  are  calculated  between  the  points  in  the  order 
that  they  are  inputted.  Summing  the  computed  distances,  the  cumulative  total  distance  to 
the  ith  point  (for  l  <  i  <  n)  is  printed  out . 

Since  straight  lines  are  drawn  between  points  by  a  computer's  plotter,  two  end  points 
are  insufficient  to  represent  a  great -circle  distance.  Therefore,  each  distance  is  broken 
down  into  shorter  line  segments,  60  nautical  miles  (1  )  in  length.  If  the  distance  is  not 
evenly  divisible  by  60,  the  last  or  end  segment  will  be  less  than  60  n.mi.  The  coordinates 
of  the  points  at  each  division  are  computed  and  then  plotted  with  the  original  end  points. 

SUBROUTINE  GCLINE  is  most  commonly  used  to  determine  a  great -circle  distance 
between  two  cities.  Inputting  the  coordinates  of  the  cities,  the  distance  is  computed  in 
nautical  miles  and  plotted.  From  the  projection  plot,  the  direction  of  the  shortest  distance 
between  the  cities  can  be  visualized.  Any  number  of  great -circle  distances  may  be  illus¬ 
trated  on  the  same  plot. 

GCLINE  may  also  be  used  to  calculate  great-circle  distances  for  consecutive  points. 

If  the  points  are  the  successive  positions  of  a  ship  sailing  between  the  first  and  last  points, 
the  total  distance  traveled  would  equal  the  sum  of  great-circle  distances. 

The  following  two  examples  are  possible  applications  of  GCLINE.  Accompanying  each 
example  is  a  listing  of  the  program,  illustrating  the  sequence  of  the  call  statements,  com¬ 
puter  printouts,  and  appropriate  plot. 
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EXAMPLE  1 


Great-circle  distances  have  been  computed  and  plotted  between  the  cities  below.  A 
hemisphere  plot  was  required  to  illustrate  these  long  distances.  Centered  at  55.0  north 
and  145.0  west,  the  plot  was  drawn  without  the  customary  meridians  and  parallels. 

•  Washington,  D.C.  to  Moscow,  USSR 

•  Los  Angeles,  California,  to  Peking,  China 
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PROGRAM  EX  AMPLE  1  1 

01  ME NS I  ON  N0PTS<25>  ,XLAT(2> ,XL0N<2>  2 

3 

4 

C  A  CALL  IS  FIRST  MAOE  TO  THE  SUBROUTINE  GLOBEMAP  TO  GENERATE  A  FULL  5 

C  GLOBE  PROJECTION.  TO  SHOW  BOTH  GREAT-CIRCLE  OISTANCES  .  THE  PLOT  6 

C  IS  CENTERED  AT  THE  POINT  (55.0  NORTH  ,  145.0  WEST).  7 

a 

9 

CALL  GLOBEMAP  (  55 . 0 , -145 . 0 .  0 .  0 , NO PT S ,  10 

*  41HWASH.  O.C.  TO  MOSCOW  ANO  L • A .  TO  PEKING, 41,  11 

*  .31HCENTERED  TO  SHOW  BOTH  01  S T ANCE S,  3 1 ,  12 

*  18HPL0T  FOR  EXAMPLE  1,18,1,6.0,6.0)  13 

14 

15 

C  STORE  THE  COORDINATES  OF  WASHINGTON  ,  D.C.  (33.55  NORTH  ,  77.0  WEST)  16 

C  ANO  MOSCOW  ,  USSR  (55.45  NORTH  ,  37.42  EAST).  THEN  CALL  GCLINE  TO  17 

C  COMPUTE  THE  GREAT-CIRCLE  OISTANCE.  18 

19 

20 

XLAT ( 1) =38. 55  21 

XL0N(1) =-77.00  22 

XLAT(2I =55.45  23 

XL0N(2) =37.42  24 

25 

CALL  GCLINE  ( XL  A T , XLON, 2 , NO PTS)  26 

27 

28 

,C  REPEAT  CALCULATION  FOR  LOS  ANGELES  (34.06  NORTH  ,  118.25  WEST)  TO  29 

C  PEKING  ,  CHINA  (39.55  NORTH  ,  116.25  EAST).  30 

31 

32 

XLAT(l) =34.  06  33 

XL0N(1)=-118.25  34 

XLAT(2) =39. 55  35 

XL0N(2I =116.25  36 

37 

CALL  GCLINE  ( XL  A T , XLON , 2 , NOPTS )  38 

39 

40 

STOP  41 

ENO  42 
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GREAT  -  CIRCLE  DISTANCE  EOUALS  4.244.16  <N.M.» 
38.55  NORTH  77.00  WEST  TD  55.45  NORTH  37.42  EAST 


LATITUDES  AND  LONGITUDES 
OF  DISTANCE  IN  ONE-DEGREE  INCREMENTS 


NO.  LAT 

• 

LONG. 

NO.  LAT, 

LONG 

• 

ND.  LAT 

• 

LONG 

• 

1-  38.6 

N, 

77.0 

W 

2-  39.4 

N» 

76.3 

W 

3-  40.2 

N, 

75.6 

W 

4-  41.0 

N, 

74.8 

w 

5-  41.9 

N, 

74.1 

M 

6-  42.7 

N, 

73.3 

w 

7-  43.5 

N, 

72.5 

H 

8-  44.3 

N* 

71.7 

H 

9-  45.1 

N, 

70.8 

w 

10-  45.9 

N, 

69.9 

W 

11-  46.7 

N, 

69.0 

W 

12-  47.4 

N, 

68.1 

w 

13-  48.2 

N, 

67.2 

w 

14-  49.0 

N  « 

66.2 

w 

15-  49.7 

N, 

65.2 

H 

16-  50.5 

N. 

64.1 

w 

17-  51.2 

N, 

63.1 

H 

18-  51.9 

N, 

62.0 

H 

19-  52.7 

N, 

60.8 

w 

20-  53.4 

N, 

59.6 

H 

21-  54.3 

N, 

58.4 

W 

22-  54.7 

N, 

57.2 

w 

23-  55.4 

N, 

55.9 

W 

24-  56.0 

N, 

54.5 

W 

25-  56.7 

N, 

53.1 

w 

26-  57.3 

N, 

51.7 

H 

27-  57.9 

N, 

50.2 

H 

28-  58.5 

N, 

48.6 

w 

29-  59.0 

N, 

47.3 

H 

30-  59.6 

N, 

45.4 

W 

31-  60.1 

N, 

43. 7 

w 

32-  60.6 

N, 

42.0 

H 

33-  61.1 

N, 

40.1 

W 

34-  61.5 

N, 

38.3 

w 

35-  62.0 

N, 

36.4 

H 

36-  62.4 

N, 

34.4 

H 

37-  62. 7 

N, 

32.4 

w 

38-  63.1 

N, 

30.3 

H 

39-  63.4 

N, 

28.2 

W 

40-  63. 7 

N, 

26.1 

w 

41-  63.9 

N, 

23.9 

w 

42-  64.1 

N, 

21.6 

H 

43 -  64 • 3 

N, 

19.4 

w 

44-  64.5 

N, 

17.1 

W 

45-  64.6 

N, 

14.8 

W 

46-  64.6 

N, 

12.5 

H 

47-  64.7 

N, 

10.1 

W 

48-  64.7 

N, 

7.8 

W 

49-  64.6 

N, 

5.4' 

H 

5C -  64.5 

N, 

3.1 

w 

51-  64.4 

N, 

0.8 

H 

52-  64.3 

N, 

1.5 

E 

53-  64.1 

N . 

3  .7 

E 

54-  63.9 

N, 

6.0 

E 

55-  63.6 

N, 

8.1 

E 

56-  63. 3 

N, 

10.3 

C 

57-  63.0 

N, 

12.4 

E 

58-  62. 7 

N, 

14.4 

E 

59-  62.3 

N, 

16.5 

E 

60-  61.9 

N, 

18.4 

E 

61-  61.5 

N, 

20.3 

E 

62-  61.0 

N, 

22.2 

E 

63-  60.5 

N, 

24.0 

E 

64-  60.0 

N, 

25.7 

E 

65-  59.5 

N, 

27.4 

E 

66-  59.0 

N, 

29.0 

E 

67-  58.4 

N, 

30.6 

E 

68-  57.8 

N, 

32.1 

E 

69-  57.2 

N, 

33.6 

E 

70-  56.6 

N, 

35.0 

E 

71-  55.9 

N, 

36.4 

E 

72-  55.5 

N, 

37.4 

E 
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GREAT  -  CIRCLE  OISTANCE  EQUALS  5,449.26  IN.M.) 
34.06  NORTH  118.25  WEST  TO  39.55  NORTH  116.25  EAST 


L ATT  TUOES  AND  LONGITUOES 
OF  OISTANCE  IN  ONE-DEGREE  INCREMENTS 


NO.  LAT. 

LONG. 

NO.  LAT. 

LONG. 

NO.  LAT 

• 

LONG. 

1-  34.1 

N, 

118.2 

w 

2-  34.8 

N, 

119.3 

w 

3-  35.6 

N, 

119.  8 

W 

4-  36.4 

N, 

120.6 

w 

5-  37.1 

N, 

121.4 

W 

6-  37.9 

N, 

122.2 

w 

7-  38.6 

N, 

123.1 

w 

8-  39.4 

N, 

123.9 

W 

9-  40.1 

N, 

124.  3 

w 

10-  40.8 

N, 

125.7 

w 

11-  41.6 

N, 

126.6 

w 

12-  42.3 

N, 

127.6 

w 

13-  43.0 

N, 

128.5 

w 

14-  43.7 

N, 

129.5 

w 

15-  44.4 

N, 

130.5 

w 

16-  45.1 

N, 

131.5 

w 

17-  45.7 

N, 

132.6 

w 

18-  46.4 

N, 

133.7 

w 

19-  47.0 

N, 

134.8 

w 

20-  47.7 

N, 

135  .9 

w 

21-  48.3 

N, 

137.1 

w 

22-  48.9 

N, 

138.3 

w 

23-  49.5 

N, 

139.5 

w 

24-  50.1 

N, 

140.7 

w 

25-  50.7 

N, 

142. C 

w 

26-  51.3 

N, 

143.3 

w 

27-  51.8 

N, 

144.7 

w 

28-  52.4 

N, 

146.1 

w 

29-  52.9 

N, 

147.5 

w 

30-  53.4 

N, 

148.9 

w 

31-  53.8 

N, 

150 . 4 

w 

32-  54.3 

N, 

151.9 

w 

33-  54.8 

N, 

153.5 

w 

34-  55.2 

N, 

155.0 

w 

35-  55.6 

N, 

156.6 

w 

36-  56.0 

N, 

158.3 

w 

37-  56.3 

N, 

160  . 

w 

38-  56.7 

N, 

161.7 

w 

39-  57.0 

N, 

163.4 

w 

40-  57. 3 

N, 

165.2 

w 

41-  57.5 

N, 

167.0 

w 

42-  57.8 

N, 

168.8 

w 

43-  58.0 

N, 

1 7C  .6 

w 

44-  58.1 

N, 

172.5 

w 

45-  58.3 

N, 

174.  3 

w 

46-  58.4 

N, 

176.2 

w 

47-  58.5 

N, 

178.1 

w 

48-  58.6 

N, 

179.9 

F 

49-  58.6 

N, 

178.0 

E 

50-  58.7 

N, 

176.1 

E 

51-  58.6 

N, 

174.2 

f 

52-  58.6 

N, 

172.3 

E 

53-  58.5 

N  , 

170.4 

E 

54-  58.4 

N, 

168.4 

F 

55-  58.3 

N. 

166.6 

E 

56-  58.2 

N, 

164.7 

E 

57-  58.0 

N, 

162.8 

E 

58-  57.3 

N, 

161.0 

E 

59-  57.5 

N. 

159.2 

E 

60-  57.3 

N. 

157.4 

c 

61-  57.0 

N, 

155.6 

E 

62-  56.7 

N, 

153.9 

E 

63-  56.3 

N, 

152.2 

F 

64-  56.0 

N, 

150.5 

F 

65-  55.6 

N, 

148.8 

E 

66-  55.2 

N  , 

147.2 

£ 

67-  54.8 

N, 

145.7 

E 

68-  54.3 

N, 

144 .1 

E 

69-  53.9 

N, 

142.6 

E 

70-  53.4 

N, 

141.1 

r 

71-  52.9 

N, 

139.7 

E 

72-  52.4 

N, 

138.3 

E 

73-  51.8 

N, 

136.9 

E 

74-  51.3 

N, 

135.5 

E 

75-  50.7 

N  • 

134.2 

c 

76-  50.2 

N, 

132.9 

E 

77-  49.6 

N. 

131.7 

E 

78-  49.0 

N, 

130.5 

F 

79-  43. 3 

N, 

129.  3 

E 

8Q-  47.7 

N, 

128.1 

E 

81-  47.1 

N, 

127.0 

E 

82-  46.4 

N, 

125.9 

E 

83-  45.8 

N, 

124  .8 

E 

84-  45.1 

N, 

123. 7 

E 

85*  44*4 

N, 

122.7 

E 

86-  43.7 

N, 

121.7 

E 

87-  43.3 

N, 

123 .7 

E 

88-  42.3 

N, 

119.7 

E 

89-  41.6 

N, 

118.8 

E 

90-  40.9 

N, 

117.9 

E 

91-  40.2 

N, 

117.0 

E 

92-  39.6 

N, 

116.3 

E 
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EXAMPLE  2 


Listed  below  are  the  coordinates  of  the  positions  of  a  ship’s  course  through  the 
Mediterranean  Sea: 

1.  36.0  north,  7.0  west  --  commence  transit 

2.  35.9  north,  2.5  west  --  1st  course  change 

3.  38.0  north,  3.0  east  --  2nd  course  change 

4.  38.0  north,  10.5  east  --  3rd  course  change 

5.  33.5  north,  16.0  east  --  4th  course  change 

6.  34.0  north,  28.0  east  -- transit  complete. 

The  total  distance  traveled  is  determined  by  calling  SUBROUTINE  GCLINE  with  the 
six  positions. 
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PROGRAM  EXAMPLE2  1 

OIMENSI ON  NOPTS<25» ,XLAT<6> ,XLON<6»  2 

3 

4 

C  STORE  THE  SIX  POSITIONS  OF  SHIP....  5 

6 

7 

DATA  (XLAT=36.0  .35.9.38. 0.38. C«  33.5.34. 2)  S 

OATA  <XLON=-7.C , -2.5,  3 . 0 , 1C . 5, 16.0 , 26. 0 >  9 

10 

11 

C  CALL  BLOWUP  FOR  ENLARGED  PROJECTION  OF  MEDITERRANEAN  SEA  ,  CENTERED  12 

C  AT  38. 5N  ,  15. CE.  THE  BLOW-UP  RANGE  ,  MEASURED  HORIZONALLY  FROM  13 

C  THE  CENTER  OF  THE  PLOT  ,  EQUALS  1150  NAUTICAL  MILES....  14 

15 

16 

CALL  BLOWUP  (35.5,15.0,0.0, 1150. C.NOPTS,  17 

*  44HENL  ARGE  O  PROJECTION  OF  THE  MEDITERRANEAN  SEA, 44,  18 

*  20HSHOWS  COURSE  OF  SHIP, 20,  19 

*  18HPLOT  FOR  EXAMPLE  2,18.0,  9. 0,6.0  20 

21 

22 

C  CALL  GCLINE  TO  DETERMINE  THE  TOTAL  OISTANCE  TRAVELED  BY  THE  SHIP....  23 

24 

25 

CALL  GCLINE  ( XL  AT , XL  ON , 6 , NOPTS »  26 

27 

STOP  28 

ENO  29 
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PLOT  FOR  EXRMPLE  2 


dlMS  JO  3SUD03  SMOHS 
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*  GREAT  -  CIRCLE  DISTANCES  * 


01 ST  ANCE 

FROM  POINT  1 

TO 

POINT 

2  EQUALS 

218.64 

(N. 

M.  > 

NO.  LAT. 

LONG. 

NO. 

LAT. 

LONG. 

NO. 

LAT 

• 

LONG. 

1“  36*0  N « 

7.  ; 

W 

2- 

36.0  N, 

5  .8  W 

3- 

36.0 

N, 

4.  5 

W 

4-  35. 9  N, 

3.3 

w 

5- 

35.9  N, 

2.5  W 

*•* 

«  * 

*  *  * 

*# 

CUMULATIVE 

DISTANCE 

TO 

POINT 

3  EQUALS 

510.84 

(N. 

M.  ) 

NO.  LAT. 

LONG. 

NO. 

LAT. 

LONG. 

NO. 

LAT 

• 

LONG. 

1-  35.9  N, 

2.5 

W 

2- 

36.4  N, 

1.4  W 

3- 

36.8 

N, 

0.3 

w 

4-  37.2  N, 

0.9 

E 

5- 

37.6  N, 

2.0  E 

6- 

38.0 

N, 

3.0 

E 

♦  * 

* 

*  ¥*  . 

*  » 

* 

CUMULATIVE 

0 1 S  T  ANCE 

TO 

POINT 

4  EQUALS 

865.35 

(N. 

M.  > 

NO.  LAT. 

LONG. 

NO. 

LAT. 

LONG. 

NO. 

LAT 

• 

LONG. 

1-  38.0  N, 

3.0 

E 

2- 

38.0  N, 

4.3  E 

3- 

38.1 

N, 

5.5 

E 

4-  38.1  N, 

6.8 

E 

5- 

3  8*1  N « 

8.1  E 

6- 

38.0 

N, 

9.3 

E 

7-  38.0  N, 

10  .5 

E 

•  » 

*  •*  » 

»  *  • 

• 

CUMULATIVE 

DISTANCE 

TO 

POINT 

5  EQUALS 

1,245.50 

(N. 

M.  ) 

* 

NO.  LAT. 

LONG. 

NO. 

LAT. 

LONG. 

NO. 

LAT 

• 

LONG. 

1*  38*0  N  y 

10.5 

2- 

37.3  N, 

11  .4  E 

3- 

36.6 

N, 

12.3 

E 

4-  35.9  N, 

13.2 

E 

5- 

35.2  N, 

14.0  E 

6- 

34.5 

N, 

14.9 

E 

7-  33.7  N, 

15.7 

E 

8- 

33.5  N, 

16.0  E 

♦  ♦ 

* 

CUMULATIVE 

01  STANCE 

TO 

i  POINT 

6  EQUALS 

1,844.57 

(N. 

M.  ) 

NO.  LAT. 

LONG. 

NO. 

LAT. 

LONG. 

NO. 

LAT 

• 

LONG. 

1-  33.5  N, 

16.0 

E 

2- 

33.6  N, 

17.2  E 

3- 

33 . 7 

N, 

18.4 

E 

4-  33.8  N, 

19.6 

E 

5- 

33.8  N. 

20.8  E 

6- 

33.9 

N, 

22.  0 

E 

7-  33.9  N, 

23.2 

E 

8- 

34.0  N, 

24.4  E 

9- 

34.0 

N, 

25.6 

E 

10-  34.0  N, 

26.8 

E 

11- 

34.0  N, 

28.0  E 
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RANGE  CIRCLES 


A  range  circle  is  computed  and  drawn  by  inputting  the  coordinates  of  its  center  and 
its  radius.  The  circle  is  divided  into  360  one-degree  arcs,  and  the  coordinates  at  each 
division  are  computed  beginning  at  true  north  clockwise  through  360  degrees.  The  coor¬ 
dinates,  which  are  used  to  draw  the  circle  on  a  projection  map,  are  also  printed  out  in 
one-degree  or  five-degree  increments.  The  computer  printouts  in  this  section  display 
the  coordinates  every  five  degrees.  The  surface  area  of  the  range  circle  is  also  com¬ 
puted  and  displayed. 

A  range  circle  is  referred  to  as  a  primary  or  secondary  range.  Executing  a  call  to 
either  of  the  subroutines,  GLOBEMAP  or  BLOWUP,  a  primary  range  is  generated  about  the 
center  point  of  the  projection.  Limited  to  one  per  plot,  a  primary  range  is  used  in  surface 
area  calculations  with  the  subroutines,  SECRANGE  and  SURFAREA. 

SUBROUTINE  SECRANGE  is  called  for  additional  secondary  ranges  radiating  from 
any  point  on  the  earth's  surface.  A  secondary  range  is  tested  to  determine  whether  it 
intersects  the  primary  range.  If  so,  the  following  are  computed: 

•  intercept  points; 

•  the  area  of  intersection; 

•  complement  surface  areas; 

•  total  surface  area  covered  by  both  ranges. 

The  primary  and  secondary  ranges  have  upper  limits  of  5,  400  and  10,  800  nautical 
miles,  respectively.  A  primary  range  is  set  equal  to  zero  if  not  desired  when  calling  the 
subroutines,  GLOBEMAP  and  BLOWUP. 

The  following  two  examples  illustrate  possible  applications  of  intersecting  range 
circles.  The  first,  example  3,  is  a  simple  two-range  intersection,  describing  an  area 
of  operation  in  the  North  Atlantic.  Example  4  introduces  the  technique  for  determining 
the  surface  area  of  multiple  range  intersections . 
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EXAMPLE  3 


The  surface  area  of  the  area  of  operation  that  is  everywhere  less  than  or  equal  to 
2,  000  n.mi.  from  New  York  City,  New  York  and  1,  200  n.mi.  of  London,  England,  is 
equal  to  the  area  of  intersection  of  the  range  circles  radiating  from  the  cities.  The  points 
of  the  resulting  area's  perimeter  have  been  outlined  on  the  following  computer  printouts. 
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PROGRAM  EXAMPLE  1 

DIMENSION  NOPT  S ( 25 )  2 

3 

4' 

C  CALL  SUBROUTINE  BLOWUP  TO  ORAM  AN  ENLARGED  PROJECTION  CENTERED  AT  5 

C  NEW  TORK  (40. TIN  ,  74. CW)  ANO  A  RANGE  CIRCLE  EDUAL  TO  2000  N.  MILES.  6 

7 

a 

CALL  BLOWUP  (40 . 71, -74. 0, 2000 .0 , 240 0 .0 , NOPTS,  9 

*  37HCENTERED  AT  NEW  YORK  CITY  ,  MEW  YORK  ,37,  10 

*  41HRANGES  EOUAL  1200  AND  2000  NAUTICAL  MILES, 41,  11 

*  18HPL0T  FOR  EXAMPLE  3,18,  12 

*  1,6.0, 6.0)  13 

14 

15 

C  DRAW  A  RANGE  CIRCLE  EOUAL  TO  1200  N.  MILES  ABOUT  LONDON  ,  ENGLAND  16 

C  (51. -.BN  ,  1.00W)  .  17 

18 

19 

CALL  SECPANGE  (51 .49,-1 . 0 , 1200. 0, NOPTS)  20 

21 

STOP  22 

END  23 
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PLOT  FOR  EXAMPLE  3 


CENTERED  AT  NEW  YORK  CITY  >  NEW  YORK 


FIG.  3:  TWO  INTERSECTING  RANGE  CIRCLES 
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PRIMARY  RANGE  EQUALS  2,01)0.00  NAUTICAL  MILES 
CENTEREO  AT  40.71  NORTH,  74.00  WEST 


SURFACE  AREA  OF  ZONE  EQUALS  1 2 . 2.T  2 , 42  3 . 4?  SO.  NAUTICAL  MILES 


LATTTUOES  AND  LONGITUDES 

DESCRIBING  TRACK  OF  RANGE  IN  INCREMENTS  OF  FIVE  DEGREES 
(BEGINNING  AT  TRUE  NORTH  AND  RGTATEQ  CLOCKWISE  THR.OUGH  360  DEGREES) 


EGREE 

LATITUDE 

LONGI TUOE 

DEGREE 

LATITUDE 

LONGITUDE 

0- 

74.04 

NORTH 

74.00 

WEST 

5- 

73 .72 

NORTH 

64.17 

WEST 

1  35- 

7.47 

NORTH 

76.  77 

WEST 

10- 

72.77 

NORTH 

55.20 

WEST 

190- 

7. 74 

NORTH 

79.53 

WEST 

15- 

71.31 

NORTH 

47.65 

WEST 

195- 

8.20 

NORTH 

82.26 

WEST 

20- 

69.45 

NORTH 

41.53 

WEST 

2  00- 

3.  83 

NORTH 

84.96 

WEST 

25- 

67.29 

NORTH 

3  7.03 

WEST 

205- 

9.64 

NORTH 

87.62 

WEST 

30- 

64.91 

north 

33.61 

WEST 

210- 

10.61 

NORTH 

90.23 

WEST 

3  5- 

62.39 

NORT  H 

31.15 

WEST 

215- 

11.75 

NORTH 

92.  73 

WEST 

40- 

59.7  7 

NORTH 

29.45 

WEST 

220- 

13.05 

NORTH 

95.26 

WEST 

45- 

57.08 

N  0  °  T  H 

28.35 

west 

225- 

1-.  50 

NORTH 

97.66 

WEST 

50- 

54.36 

NORTH 

27.75 

WEST 

2  3  0- 

16.09 

NORTH 

99.  98 

WEST 

55- 

51.61 

NORTH 

27.54 

WEST 

235- 

IT.  82 

NO°TH 

102.22 

WEST 

60- 

48.87 

NORTH 

27.66 

WEST 

240- 

19.67 

NORTH 

104.36 

WEST 

65- 

46.13 

NORTH 

28.0  5 

WEST 

245- 

21.65 

NORTH 

106.40 

WEST 

70- 

43.42 

NORTH 

28.69 

WEST 

250- 

23.73 

NORTH 

lOd. 34 

WEST 

75-  • 

40.75 

NORTH 

29.52 

WEST 

255- 

25.92 

NORTH 

110.17 

WEST 

80- 

33.12 

NORTH 

30.54 

WEST 

260- 

23.2  0 

NORTH 

111.33 

WEST 

85- 

35 .54 

NORTH 

31.72 

WEST 

265- 

30.57 

north 

113.43 

WEST 

90- 

33.0  2 

NORTH 

33.05 

WEST 

270- 

33.02 

NORTH 

114.95 

WEST 

95- 

30.57 

NORTH 

34.52 

WEST 

275- 

35.54 

NORTH 

116.28 

WEST 

100- 

28.20 

NORTH 

36.12 

WEST 

2  80- 

33.12 

NO=TH 

117.46 

WEST 

105- 

25.92 

NORTH 

37.83 

WEST 

285- 

40. 75 

NORTH 

118.48 

WEST 

110- 

23 . 73 

NORTH 

39.66 

WEST 

290- 

43.4? 

NORTH 

119.31 

WEST 

115- 

21.65 

NORTH 

41.60 

WEST 

295- 

46.13 

NORTH 

119.95 

WEST 

120- 

19.67 

•NORTH 

43.64 

WEST 

300- 

43.87 

NORTH 

120.34 

WEST 

125- 

17 .82 

NORTH 

-5.78 

WEST 

305- 

51.61 

NO=>TH 

120 . 46 

WEST 

130- 

16.09 

NORTH 

48.02 

WEST 

3  10- 

54.  36 

NORTH 

120. '25 

WEST 

13  5- 

14.50 

NORTH 

50.34 

WEST 

315- 

57.03 

NORTH 

119.65 

WEST 

140- 

13.05 

NORT  H 

52.74 

WEST 

320- 

59.77 

NQD  T  H 

11.3.55 

WEST 

145- 

11.75 

NORTH 

55.22 

WEST 

32  5- 

62.39 

NORTH 

116.85 

WEST 

150- 

10.61 

NORTH 

57.  77 

WEST 

3  30- 

64.91 

NORTH 

114.39 

WEST 

155- 

9.64 

NORTH 

60.38 

WEST 

335- 

67.29 

NORTH 

110.97 

WEST 

160- 

8.83 

NORTH 

63.04 

WEST 

34  0- 

6  9.45 

NORTH 

106. 37 

WEST 

16  5- 

8.20 

NORTH 

65.74 

WEST 

345- 

71.31 

NORTH 

100.35 

WEST 

170- 

7.74 

NORTH 

68.47 

WEST 

350  - 

72.77 

NORTH 

92.80 

WEST 

175- 

7.47 

NORTH 

71.23 

WEST 

355- 

73.  72 

NORTH 

33. 33 

WEST 

180- 

7.38 

NORTH 

74.00 

WEST 

360- 

74.04 

NORTH 

74.00 

WEST 
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SECONDARY  RANGE  EQUALS  1,200.00.  NAUTICAL  MILES 
CENTtREO  AT  51. 4<J  NORTH,  1.00  WEST 


SURFACE  AREA  OF  ZONE  EQUALS  4,454,199.06  SO.  NAUTICAL  MILES 


LATITUOES  ANO  LONGITUDES 

DESCRIBING  TRACK  Oc  RANGE  IN  INCREMENTS  OF  FIVE  DEGREES 
(BEGINNING  AT  TRUE  NORTH  ANO  ROTATED  CLOCKWISE  THROUGH  360  DEGREES) 


OEGREE 

LATITUDE 

LONG  I TUOE  OEGREE 

LATITUDE 

LONGITUDE 

0- 

71.49 

NORTH 

i .  a  a 

WEST 

5- 

71.34 

NORTH 

4.35 

EAST 

185- 

31.54 

NORTH 

3.00 

WEST 

10- 

70.91 

NORTH 

9.47 

EAST 

190- 

31.71 

NORTH 

5.00 

WEST 

15- 

70.22 

NORTH 

14.17 

FAST 

195- 

31.98 

NORTH 

6.99 

WEST 

20- 

69 . 30 

NORTH 

13.32 

EAST 

200- 

32.36 

NORTH 

8.96 

WEST 

25- 

65.17 

NORTH 

21.38 

EAST 

205- 

32.84 

NORTH 

10.91 

WEST 

30- 

66.89 

NORTH 

24.83 

EAST 

2  10- 

33.4  3 

NORTH 

12.82 

WEST 

35- 

65.47 

NORTH 

27.20 

EAST 

215- 

34.12 

NORTH 

14.71 

WEST 

40- 

63.95 

NORTH 

29.05 

EAST 

220- 

34.  90 

NORTH 

16.55 

WEST 

45- 

62.36 

NORTH 

3  0.-2 

EAST 

225- 

35.78 

NORTH 

18.34 

WEST 

50- 

60.71 

NORTH 

31.35 

EAST 

230- 

36.76 

NORTH 

20.09 

WEST 

55- 

59.03 

NORTH 

31 .99 

EAST 

235- 

37.82 

NORTH 

21.77 

WEST 

60- 

57  .33 

NORTH 

32.28 

EAST 

240- 

38.96 

NORTH 

23.39 

WEST 

65- 

55.62 

north 

32.29 

EAST 

2  4  3- 

40.19 

NORTH 

24.94 

WEST 

70- 

53.92 

NODT  H 

32.07 

EAST 

250- 

41 . 49 

NORTH 

26.41 

WEST 

75- 

52.23 

NORTH 

31.04 

EAST 

255- 

42.86 

NORTH 

27.79 

WEST 

80- 

50.56 

NORTH 

31.02 

EAST 

260- 

44.29 

NORTH 

29.0  7 

west] 

85- 

48.93 

NORTH 

30.24 

EAST 

265- 

45.79 

NORTH 

30.25 

west] 

9  0- 

47.33 

NORTH 

2  9.31 

EAST 

270- 

47.33 

NORTH 

31.31 

WEST 

95- 

45.79 

NORT  H 

23.25 

EAST 

275- 

48.93 

NORTH 

32.24 

WEST] 

100- 

44.29 

NORTH 

27.07 

EAST 

280- 

50.56 

NORTH 

33.0  2 

WEST 

105- 

42,66 

NORTH 

25.79 

EAST 

285  - 

52.23 

NORTH 

33.64 

WEST 

110- 

41.49 

NORTH 

24.41 

EAST 

290- 

53.92 

NORTH 

34.07 

WEST 

115- 

40.19 

NORTH 

22.  94 

EAST 

295- 

55.62 

NORTH 

34.29 

WEST 

120- 

38.96 

NORTH 

21.39 

EAST 

300- 

57.  33 

NO=TH 

34.28 

WEST 

125- 

37.82 

NORTH 

19.77 

EAST 

305- 

59.0  3 

NORTH 

33.99 

WEST 

130- 

36.76 

NORTH 

13.0  9 

EAST 

3  10- 

60.71 

NORTH 

3  3.39 

W  E  3  T 

135- 

35.73 

NC  =  TH 

16.34 

EAST 

315- 

62.36 

NO=>TH 

32.  *2 

WEST 

140- 

3  -♦ .  9  0 

NORTH 

14.55 

EAST 

320- 

63.95 

NORTH 

31.05 

WEST 

145- 

34.12 

NORTH 

12.71 

EAST 

325- 

65.47 

NORTH 

29.20 

WEST 

150- 

33.43 

NORTH 

10  .  82 

EAST 

3  30  - 

66.89 

NORT  H 

26.33 

WEST 

155- 

32.84 

NORTH 

8.91 

EAST 

335- 

68.17 

NORTH 

23.83 

WEST 

160- 

32.36 

NORT  H 

6.  96 

EAST 

3  40- 

69.30 

NORTH 

20.32 

WEST 

165- 

31.93 

NORTH 

4 . 99 

EAST 

345- 

70.22 

NORTH 

16.17 

WEST 

170- 

31.71 

NORTH 

3.00 

EAST 

350- 

70.91 

NORTH 

11.47 

WEST 

175- 

31.54 

NORTH 

1.00 

EAST 

355- 

71.34 

NORTH 

6.35 

WEST 

180- 

31.49 

NORTH 

1.00 

WEST 

360- 

71.49 

NO°T  h 

1.00 

WEST 

SECON 0 AR Y  PANG-  INTERSECTS  PRIMAPY  RANGE  AT 
(6  3.  to  NORT*.  H  .  7 9  -v-rm  ANO  143.75  NO->TH.  2°,  60  WEST] 

SURFACE  AREA  GF  INTERSECTION  EQUALS  177,381. 7h  SO.  NAUTICAL  MILES 


SURFACE 

AREA 

OF 

REMA  INING 

PRIMARY 

ZONE 

EQUALS 

12,055,041.75 

SQ. 

NAUTICAL 

MILES 

SURFACE 

AREA 

OF 

REMAINING 

SECONDARY 

ZONE 

EQUALS 

4,306,817.33 

SO. 

NAUTICAL 

MILES 

TOTAL 

COVE 

RAGE  AREA 

OF  BOTH  RANGES 

EQUALS 

16,539,240. 8? 

SQ. 

NAUTICAL 

MILES 
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EXAMPLE  4 


For  three  range  circles  centered  at  Canton,  Peking,  and  Tokyo,  and  having  radii  equal 
to  750,  800  and  900  nautical  miles,  respectively,  find  the  surface  areas  of  double  and  triple 
coverages. 

From  the  computer  printouts  the  three  areas  of  intersection  equal  453,  293;  25,  881; 
and  492,  861  n.mi.l  Subtracting  the  area  of  triple  coverage,  2  which  equals  17,  995  n.mi., 
from  each  of  the  above  areas,  the  surface  area  of  exactly  double  coverage  equals  918,  050 
n.mi.  The  perimeter  points  of  the  triple  intersection  have  been  outlined  on  the  printouts. 


1This  problem  required  an  additional  plot  to  determine  the  area  of  intersection  between 
Peking  and  Tokyo  circles.  Only  the  plot  centered  at  Canton  is  presented. 

2 

See  example  6,  page  37,  for  the  surface  area  solution  of  the  triple  coverage. 
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PROGRAM  EXAMPLE* 
DIMENSION  NOP TS ( 25 1 


C  CALL  SU8R0UTI NE  BLOWUP  TO  DRAW  AN  ENLARGEO  PROJECTION  CENTEREO  AT 
C  CANTON  ,  CHINA  (23.03N  ,  113. 20EJ.  PRIMARY  RANGE  EQUALS  750  NAUTICAL 
C  MILES . 


CALL  BLOWUP  ( 23 . 03, 113. 20, 750 .0 .2200 .0 , NOPTS, 

39 HRANGES  EQUAL  750  ,  80  0  AND  900  N.  MILES, 39, 
45HFR0M  CANTON  ,  PEKING  AND  TOKYO  ,  RESPEC T I VEL Y , 45 , 
18HPLOT  FOR  EXAMPLE  4,18, 

1,6. 0,5.0) 


C  DRAW  RANGE  CIRCLE  OF  800  N. MILES  FROM  PEKING  ,  CHINA.... 


CALL  SECRANGE  ( 3 9. 5 5 , 1 1 6 . 25 , 80 0 . 0 , NOPTS ) 


C  ORAW  RANGE  CIRCLE  OF  900  N. MILES  FPOM  TOKYO  ,  JAPAN.... 


CALL  SECRANGE  ( 35 . 6 7 . 13 9 . 75 , 9 0 Q . 0 , N OP TS ) 

STOP 

ENO 


1 

2 

3 

•4 

5 

6 
7 

a 

9 

10 

11 

12 

13 

14 

15 

16 
17 
19 

19 

20 
21 
22 

23 

24 

25 

26 
27 
29 
29 
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FROM  CANTON  >  PEKING  AND  TOKYO  >  RESPECTIVELY 


RANGES  EQUAL  750  >  800  AND  900  N.  MILES 
FIG.  4:  THREE  INTERSECTING  RANGE  CIRCLES 
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.PRIMAPY  RANGE  EQUALS  750.00  NAUTICAL  MILES 
CENTEREO  AT  23.08  NORTH,  113.20  EAST 


SURFACE  AREA  OF  ZONE  EQUALS  1,76?, 527. 64  SQ.  NAUTICAL  MILES 


LATITUQES  ANO  LONGITUOES 

0ESCRI9ING  TRACK  OF  range  IN  INCREMENTS  of  FINE  degrees 
(0EGINNING  AT  TRUE  NORTH  AND  kOTATEO  CLOCKWISE  THROUGH  360  DEGREES) 


OEGREE  LATITUDE  LONG  I TIJOE 


0- 

35.58 

NORTH 

113.20 

EAST 

5- 

35.53 

NORTH 

114.53 

EAST 

10- 

35.37 

NORTH 

113.34 

EAST 

15- 

35.10 

NORTH 

117.13 

EAST 

20- 

34. 74 

NORTH 

113.37 

east 

25- 

34.28 

NORTH 

119.56 

EAST 

30- 

33.72 

NORTH 

12  0.68 

EAST 

35- 

33.U8 

NORTH 

121. 72 

TTsT 

40- 

32.36 

NORTH 

122.68 

EAST 

45- 

31.57 

NORTH 

123.55 

EAST 

50- 

30.71 

NORTH 

124.32 

EAST 

55- 

29.80 

NORTH 

124.99 

EA  ST 

60- 

oc 

<r 

j 

NO°  TH 

123.56 

EAST 

65- 

27.83 

NORTH 

126.0  2 

EAST 

70- 

26.80 

north 

126.37 

EAST 

75- 

25.74 

NORTH 

126.62 

EAST 

80- 

2-  .66 

NORTH 

126.77 

EAST 

85- 

23.58 

NORTH 

126.31 

EAST 

90- 

22.50 

NO°TH 

126.75 

EAST 

95- 

21.43 

NORTH 

12^.59 

EAST 

100- 

20.37 

NORTH 

126.34 

EAST 

105- 

19.34 

NOR  TH 

126.00 

EAST 

110- 

18.34 

NC-TH 

12-3.5  7 

EAST 

115- 

17.37 

NORTH 

12  5.06 

EAST 

120- 

lo  .  4  5 

NORTH 

12-.. -.7 

EAST 

125- 

15.58 

NORTH 

123.81 

EAST 

130- 

14.76 

•NORTH 

123.07 

EAST 

135- 

14.00 

NORTH 

1  22.28 

EAST 

140- 

13.31 

NORTH 

121.42 

EAST 

145- 

12.  69 

NOET  H 

120.51 

FAST 

150- 

12.14 

NORTH 

119.56 

EAST 

155- 

11.67 

NORTH 

113.56 

EAST 

160- 

11 .2d 

NORTH 

117.53 

EAST 

165- 

10.98 

NORTH 

116.  -.7 

EAST 

170- 

10.76 

NORTH 

115.39 

EAST 

175- 

10.62 

NORTH 

114.30 

EAST 

180- 

10.53 

NORTH 

113.20 

EAST 

DEGREE  LATITUDE  LONGITUDE 


185- 

10.62 

NO=>TH 

112.19 

EAST 

190- 

10.76 

NORTH 

111.01 

EAST 

195- 

10.98 

NORTH 

109.93 

EAST 

200- 

11.28 

NORTH 

108.87 

fast 

205- 

11.67 

NORTH 

107.84 

EAST 

210- 

12.14 

NORTH 

106.84 

EAST 

215- 

12.69 

NO^TH 

105.89 

EAST 

220- 

13.31 

NORTH 

104.98 

EAST 

225- 

14.0  0 

NORTH 

104.12 

EAST 

230- 

14.76 

NORTH 

103.33 

EAST 

2  35- 

15.58 

NORTH 

102.59 

EAST 

240- 

15.45 

NO^TH 

101.93 

EAST 

245- 

17.37 

NORTH 

101.34 

EAST 

250  - 

18.34 

NORTH 

100.83 

EAST 

255- 

19.34 

NORTH 

100.40 

FAST 

260- 

20.37 

NORTH 

100.06 

EAST 

265- 

21.43 

NORTH 

99.31 

EAST 

270- 

22.50 

NO3  TH 

99.65 

EAST 

275- 

23.58 

NORTH 

99.59 

EAST 

230- 

24.66 

NORTH 

99.63 

EAST 

285- 

25.74 

NO  ? T H 

99.73 

EAST 

290- 

25.30 

NORTH 

100.03 

FAST 

295- 

27.83 

NORTH 

100.33 

EAST 

3  00- 

28.83 

NORTH 

100.34 

EAST 

305- 

29.8  0 

NORTH 

101.41 

EAST 

310- 

30.71 

NORTH 

102.08 

EAST 

315- 

31.57 

NORTH 

102. 35 

EAST 

320- 

32.36 

NORTH 

103.72 

EAST 

325- 

33.08 

NORTH 

104.68 

EAST 

330- 

33.72 

NORTH 

105.72 

EAST 

3  35- 

34.28 

NORTH 

106. 34 

EAST 

340- 

34.74 

NORTH 

103. 03 

EAST 

345- 

35.10 

NORTH 

109.27 

EAST 

3  50- 

35.37 

NORTH 

110.56 

EAST 

355- 

35.53 

NORTH 

111.37 

EAST 

3  6  0  - 

35.58 

NO15  TH 

113.20 

EAST 
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SECONOARY  RANGE  EQUALS  800.00  NAUTICAL  NILES 
CENTEREO  AT  39. 5E  NORTH,  116.25  CAST 


SURFACE  AREA  OF  ZONE  EQUALS  2,304,263.26  SQ.  NAUTICAL  NILES 

»*¥**»»»*» 

LATITUDES  ANO  LONGITUDES 

0ESCRI3ING  TRACK  OF  RANGE  IN  INCREMENTS  OF  FIVE  OEGREES 
(BEGINNING  AT  TRUE  NORTH  ANO  ROTATED  CLOCKWISE  THROUGH  360  OEGREES) 


OEGREE 

LAT  IT 'JOE 

LONGITUOE 

OEGREE 

LATITUOE 

LONGITUOE 

0- 

52.88 

NCRT  H 

llo.25 

EAST 

5- 

52.82 

NORTH 

113.16 

EAST 

1  85- 

26.26 

NORTH 

114.97 

EAST 

10- 

52.63 

NORTH 

120.03 

EAST 

190  - 

26.39 

NORTH 

113.69 

EAST 

15- 

52.31 

NORTH 

121.85 

EAST 

195- 

26.60 

north 

112.42 

EAST 

20- 

51.83 

NORTH 

123.59 

EAST 

200- 

26.90 

NORTH 

111.18 

EAST 

25- 

51  .33 

NORTH 

125.22 

EAST 

205- 

27.29 

NORTH 

109.95 

EAST 

30- 

50.63 

NORTH 

126. 73 

EAST 

210- 

27.75 

NORTH 

108.76 

EAST 

35- 

49.9  3 

NORTH 

123.11 

EAST 

215- 

28.29 

NORTH 

107.61 

EAST 

40- 

49.10 

north 

129.33 

EAST 

220  - 

28.91 

NORTH 

106. 50 

EAST 

4  5- 

48.19 

MORT  H 

130  .41 

EAST 

225- 

29.59 

NORTH 

105.44 

EAST 

50- 

47.21 

north 

131.32 

EAST 

2  30- 

30.35 

NORTH 

104.44 

EAST 

55- 

46.19 

NORTH 

132.08 

EAST 

235- 

31.17 

NORTH 

103.49 

EAST 

60- 

45.11 

NORTH 

1 32.69 

EAST 

240- 

32.05 

NORTH 

102.62 

EAST 

65- 

44.0  1 

NORTH 

133.14 

EAST 

2  45  - 

32.99 

NORTH 

101.32 

EAST 

70- 

42.88 

NORTH 

133.45 

EAST 

250- 

33.97 

NORTH 

101.10 

EAST 

75- 

41.73 

NORTH 

133.62 

EAST 

255- 

35.00 

NORTH 

100.47 

EAST 

80- 

40.5  3 

NORTH 

133.65 

EAST 

260- 

36.07 

NORTH 

99.93 

EAST 

85- 

39.43 

NORTH 

133.55 

EAST 

265- 

37.16 

NORTH 

99.49 

EAST 

90- 

3  8.29 

NORTH 

133. 34 

EAST 

270- 

33.29 

NORTH 

99.16 

EAST 

9  5- 

37.16 

NORTH 

133.01 

EAST 

275- 

39.  43 

NORTH 

9  3.95 

EAST 

100- 

36.07 

NORTH 

132.57 

EAST 

2  80- 

40.53 

NORTH 

99.35 

EAST 

10  5- 

35.00 

NORTH 

132.03 

EAST 

2  85- 

41.73 

NORTH 

99. 98 

EAST 

110- 

33.97 

NORTH 

131.40 

EAST 

290- 

42.88 

NORTH 

99.05 

EAST 

115- 

32.99 

NORTH 

130.68 

EAST 

295- 

44.01 

NORTH 

99.36 

EAST 

120- 

32. C  5 

NORTH 

129. 88 

EAST 

3  00- 

45.11 

NO  =  TH 

99.  91 

EAST 

12  5- 

31.17 

NORTH 

12  9.01 

EAST 

3  0  5- 

4  6.19 

NORTH 

100.42 

EAST 

130- 

30.35 

NORTH 

128.06 

EAST 

310  - 

47.21 

NORTH 

101.18 

EAST 

135- 

29.59 

NORTH 

127.06 

EAST 

315- 

43.19 

NORTH 

10  2.09 

EAST 

140- 

23.91 

NORTH 

1 26.00 

EAST 

320- 

49.10 

NORTH 

103.17 

EAST 

1145- 

28.29 

*  r.7  ^  j  h 

124.89 

r.  A  S  T| 

.  325- 

49.93 

NORTH 

104.39 

EAST 

150- 

27.75 

NORTH 

123.74 

EAST 

3  30  - 

50.68 

NORTH 

105.77 

EAST 

155- 

27.29 

NORTH 

122  .55 

EAST 

3  3  5- 

51.33 

NORTH 

107.23 

EAST 

160- 

26  ,Q0 

NORT  H 

121.32 

EAST 

340- 

51.88 

NORTH 

108.91 

EAST 

165- 

26.60 

nc°th 

12  0.08 

EAST 

345- 

52.31 

NORTH 

110 .65 

EAST 

170- 

26.39 

NORTH 

118.31 

EAST 

3  50- 

52.63 

NORTH 

112.47 

EAST 

175- 

26.26 

NORTH 

117.53 

EAST 

355-. 

52.82 

NORTH 

114.34 

EAST 

180- 

26.22 

NORTH 

116.25 

EAST 

360- 

52.88 

NORTH 

116.25 

EAST 

SECONOARY  PANGE  INTERSECTS  PR  I  MA°  Y  RANGE  AT 
31.69  NORTH.  102. 97  EAST  AND  123.6  V  NORTH,  125.63  ~  A  S  fl 


SURFACE  AREA 

OF  INTERSECTION 

EQUALS  453, 

292.33  SQ. 

NAUTICAL  MIL 

ES 

SURFACE 

AREA 

OF  REMAINING 

PRIMARY  ZONE 

EQUALS  1, 

309,234.81 

SO. 

NAUTICAL 

MIL 

SURFACE 

AREA 

OF  REMAINING 

SECONDARY  ZONE 

EQUALS  1, 

550 ,975.43 

SO. 

NAUTICAL 

MIL 

total 

COVERAGE  AREA 

OF  BOTH  RANGES 

♦  ♦  *  *  * < 

EQUALS  3* 

313,503.07 

SO. 

NAUT ICAL 

MIL 

-29- 


SECONDARY  RANGE  EQUALS  900.00  NAUTICAL  MILES 
.  CENTEREO  AT  35.67  NORTH,  139.75  EAST 

********** 

SURFACE  AREA  OF  70NE  EQUALS  2.533,609.98  SQ.  NAUTICAL  MILES 

********** 

LATITUDES  AND  LONGITUDES 

DESCRIBING  TRACK  OF  PANGE  IN  INCREMENTS  OF  FIVE  OEGREES 
(BEGINNING  AT  TRUE  NORTH  AND  ROTATED  CLOCKWISE  THROUGH  360  DEGREES) 


OEGREE  LATITUDE  LONGITUOE 


0- 

50 .67 

NORTH 

139.75 

EAST 

5- 

50  .60 

NCRIH 

141.79 

EAST 

10- 

50  .38 

NCRTH 

143. 79 

EAST 

15- 

50.03 

NORTH 

1 45 • 7  4 

EAST 

20- 

49 .54 

NORTH 

147.59 

EAST 

25- 

48.92 

NORTH 

149.33 

EAST 

30- 

48.19 

NO’TH 

150  .94 

EAST 

35- 

47.35 

NORTH 

152.41 

EAST 

40- 

46.41 

NORTH 

153.71 

EAST 

4  5- 

45.39 

NORTH 

154.96 

EAST 

50- 

44.30 

NORTH 

155.83 

EAST 

55- 

43.15 

NORTH 

15  6 .64 

EAST 

60- 

41  .94 

NORTH 

157.29 

EAST 

65- 

40.70 

NORTH 

157.77 

EAST 

70- 

39.43 

NORTH 

153.10 

EAST 

75- 

38.15 

NORTH 

15  3.29 

EAST 

80- 

36.85 

NORTH 

153 . 32 

EAST 

85- 

35.56 

NORTH 

159.23 

EAST 

90- 

34.23 

NORTH 

159.00 

FAST 

95- 

33.02 

NORTH 

157.66 

EAST 

100- 

31.79 

NORTH 

157.20 

EAST 

105- 

30.59 

NORTH 

156.63 

EAST 

110- 

29.43 

NORTH 

155.96 

EA  ST 

115- 

23.3  2 

NO=tm 

155.20 

EAST 

120- 

27.27 

north 

154.36 

EAST 

125- 

26  .27 

north 

133.43 

FAST 

130- 

25.35 

north 

152.42 

EAST 

135- 

24.49 

ncrt  h 

151.35 

EAST 

140- 

23 . 71 

north 

150.22 

EAST 

145- 

23.02 

north 

149.03 

EAST 

150- 

22.41 

•  NORTH 

147.30 

EAST 

155- 

21.89 

NORTH 

146.52 

EAST 

160- 

21.45 

NCRTH 

145.21 

EAST 

165- 

21.11 

NORTH 

143.37 

EAST 

170- 

20.37 

NORTH 

142.51 

EAST 

175- 

20 .72 

NORTH 

141.13 

EAST 

180- 

20 .67 

NORTH 

139.75 

EAST 

OEGREE 

LATI TUDE 

LONGI TUDE 

185- 

20.72 

NORTH 

138.37 

EAST 

190- 

20 .87 

NO  =  TH 

136.99 

EAST 

195- 

21.11 

NORTH 

135.63 

FAST 

200- 

21.  h5 

NORTH 

134.29 

EAST 

205- 

21.88 

NORTH 

132.93 

EAST 

210- 

22.41 

NORTH 

131.70 

EAST 

215- 

2  3.02 

NORTH 

133.47 

EAST 

220- 

23.71 

NORTH 

129.28 

EAST 

225- 

24.49 

NORTH 

123.15 

EAST 

230- 

25.35 

NORTH 

127.03 

EAST 

235- 

26.27 

NORTH 

126.07 

FAST 

240- 

27.27 

NORTH 

125.14 

EAST 

2  4  5- 

28. 32 

NORTH 

124.30 

EAST 

250- 

29.43 

NORTH 

123.54 

EAST 

255- 

30.59 

NORTH 

122 . 97 

EAST 

260- 

31.79 

NORTH 

122.30 

EAST 

265  - 

33.02 

NORTH 

121.34 

EAST 

270- 

34.29 

NORTH  . 

121.53 

CAST 

275- 

35.56 

NORTH 

121.27 

FAST 

280- 

35.85 

NORTH 

121.18 

EAST 

285- 

38.15 

NORTH 

121.21 

EAST 

290- 

39.43 

NORTH 

121.40 

EAST 

295- 

40.70 

NORTH 

121.73 

r  A  S  T 

3  00- 

41. 94 

NORTH 

122.21 

EAST 

305- 

43.15 

NORTH 

122.36 

EAST 

310- 

44.30 

NORTH 

123.67 

EAST 

3  15- 

45.39 

NORTH 

124.64 

EAST 

3  20- 

4  6*41 

NORTH 

125.79 

EAST 

325- 

47.35 

NORTH 

127.09 

EAST 

3  30- 

43.19 

NORTH 

129.56 

EAST 

335- 

49.92 

NORTH 

130.17 

EAST 

340- 

49.  54 

NORTH 

131.91 

EAST 

345  - 

50.03 

NORTH 

133.76 

EAST 

350- 

53.39 

NORTH 

135.71 

EAST 

355- 

50.60 

NORTH 

137.71 

EAST 

360- 

50.67 

NORTH 

139.75 

EAST 

*****  ****  * 

SECONDARY  CANGE  INTERSECTS  PRIMARY  RANGE  AT 
|32.  99  N  G  R  T  -  ,  121.85  QsTJ  AND  25.75  NORTH,  126.62  EAST 


SURFACE  AREA  CF  INTERSECTION  EQUALS  25,881.13  SQ.  NAUTICAL  MILES 


SURFACE  AREA 

OF  REMAINING 

1 

PRIMARY 

ZONE 

*  *  **  * 

EQUALS 

1.736, 

646.51 

SO.  NAUTICAL 

MIL 

SURFACE  AREA 

OF  REMAINING 

SECONOAR.Y 

ZCNE 

EQUALS 

2.507, 

729.85 

SQ.  NAUTICAL 

MIL 

TOTAL 

COVERAGE  AREA 

OF  OOTH  RANGES 

EQUALS 

4.270, 

256.49 

SQ.  NAUTICAL 

MIL 
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SURFACE  AREAS 


SUBROUTINE  SURFAREA  will  compute  the  surface  area  of  the  portion  of  the  earth's 
surface  enclosed  by  a  set  of  points  describing  its  perimeter.  An  area  is  assumed  to  have 
great-circle  sides  and  is  computed  in  square  nautical  miles.  The  length  of  its  perimeter 
is  also  computed. 

Surface  Area  Within  a  Range 


SURFAREA  has  the  additional  capability  to  determine  how  much  of  an  area  is  within 
a  particular  distance  of  a  specified  point.  The  procedure  is  to  intersect  the  area  with  a 
range  circle  and  calculate  the  surface  area  lying  within  the  circle.  Since  a  range  circle’s 
circumference  isn't  a  great -circle  arc,  the  part  that  intersects  the  area  is  broken  down 
into  one  degree  great-circle  line  segments.  If  an  area  has  an  irregular  shape,  two  or 
more  intersections  may  occur  with  the  range  circle,  as  illustrated  in  example  5.  The 
various  computed  areas  must  be  added  or  subtracted  to  determine  the  total  surface  area 
within  the  circle. 

Inputs 

SURFAREA' s  input  parameters  to  the  Fortran  statement, 

CALL  SURFAREA  (XLAT,  XLON,  NP,  CLAT,  CLON,  NOPTS), 
are  described  as  follows: 

XLAT,  XLON  =  arrays  of  latitudes  and  longitudes  of  NP  number  of  points. 

NP  =  number  of  points . 

CLAT,  CLON  =  coordinates  locating  the  position  of  alphanumeric  character, 
identifying  the  area . 

NOPTS  =  inhouse  options. 

Rules  and  Limitations 


•  The  accuracy  in  calculating  a  surface  area  is  dependent  upon  the  selection 
of  the  points  outlining  the  area.  A  fraction  of  a  degree  difference  in  one  of 
the  points  can  vary  the  total  area  by  several  thousand  square  nautical  miles. 

•  The  coordinates  of  the  first  and  last  points  must  be  equal  to  close  the  area. 
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•  The  exception  to  the  above  occur  when  intersecting  an  area  with  a  range 
circle.  The  first  and  last  points  can  be  selected  to  lie  just  outside  the 
circle  and  need  not  be  equal.  When  working  with  a  large  area,  only  the 
perimeter  points  inside  and  just  outside  the  range  circle  have  to  be  inputted. 

•  The  routine  that  controls  the  plotting  of  all  points  also  checks  the  points  of 
an  area's  perimeter  to  determine  whether  they  fall  outside  the  plot  border. 

If  they  do,  these  points  are  omitted  from  the  calculation.  The  resulting  sur¬ 
face  area  is  erroneous  and  usually  negative. 

•  The  size  of  an  area  has  an  upper  limit  equal  to  a  hemisphere  of  74,  355,  728 
square  nautical  miles. 

•  The  perimeter  of  an  area  must  be  represented  by  at  least  four  points. 
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EXAMPLE  5 


The  objective  of  this  example  is  to  illustrate  four  different  surface  area  calculations. 
Given  the  perimeter  coordinates  of  the  four  areas,  the  surface  areas  within  2, 000  nautical 
miles  of  Washington,  D.C.  are  calculated.  The  coordinates  of  the  perimeters  can  be 
found  in  the  program  listing. 

The  following  results  were  observed: 

•  AREA  A  was  entirely  enclosed  within  the  2,000  n.mi.  range  circle.  Its 
surface  area  and  perimeter  were  290,  288.0  square  n.mi.  and  4,171.7  n.mi., 
respectively. 

•  AREA  B  was  intersected  by  the  range  circle  yielding  a  surface  area  equal  to 
361,  348 . 4  square  n.mi.  The  perimeter  which  includes  the  part  of  the  range 
circle  intersecting  the  area  equals  4,  617.8  n.mi. 

•  AREA  C,  due  to  its  irregular  shape,  was  intersected  in  three  places.  The 
total  surface  area  within  the  circle  equals  676,045  square  n.mi. 

[C(l)  +  C(2)  -  C(3)  ]  . 

•  Since  AREA  D  was  outside  the  range,  it  was  not  plotted  and  no  surface  area 
was  calculated. 

Since  SURFAREA  will  not  plot  pass  the  primary  range,  the  perimeters  lying  outside 
have  been  drawn  by  hand.  SUBROUTINE  GCLINE  could  have  been  called  to  plot  the  entire 
areas. 
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PROGRAM  E  XA  MPLE  5 

01  MENSI  ON  N  0°  TS  (P.5)  ,  AL  A  T  (30)  ,ALON(30),9LAT(30),0L0N(30>, 
CL A  T (30) , CLON  (  30)  . OLA T (30 > ,DLQN( 30) 


C  STORE  COORDINATES  OF  PERIMETER  POINTS  OUTLINING  THE  FOUR  AREAS.... 


OATA  (ALA T= 44, 0,42.5,39.5,3ft. 0,35. 0,30. 0,27. 0,22. 0,21. 0,21. 0,23.0. 

►  25. 5, 29. 0,31. 0,3 6. 0,40. 0,43. 0.44. 0,44.0) 

OATA  (AL0N=-E2. U, -5  7. 5, -59. 0,-73. 0,-74. 0,-79. 0, -75.0,-70.0, -59.0, 

*  -65.0,-65.0,-70.0,-74.0,-75.0,-70.0,-55.0,-61.0,-61.0, 

►  -62.0) 

OATA  ( BL A T= 22. 5, 22. 0,22. 0,20. 0,1ft. 0,15. 5, 14. 0,15. 0,1?. 0,12. 0,8.0, 

*  6.5,7.5,5.01 

OATA  (BLON= -113.. 0,-111. 0,-105. 0,-107. 0,-104. 0,-100.0, -97.0, -95.0, 

*  -91. 0,-bft.n, -rt5. 0,-ftl.0, -79. 0, -78.5) 

OATA  (CLAT=58.0, 53.0,5fl. 0, 58. 0,55. 0,55. 0,55. 0,55. 0,51. 0,51. 0,51.0. 

►  51. 0,51. 0.50. 0,49. 0,47. 0,44. 0,4 0.0, 33. 0,3 2. 0,35. 0,37.0, 

»  39. 0,40. 0.41. 5, 43. 0,45. 5, 50. 0,53.0) 

DATA  (CLON= -2 1.0, -30. 0,-40. 0,-4 6.0, -46.0,-40.0, -30.0, -29.0, -29.0, 

*  -30.0,-40.0,-50.0,-52.0,-51.0.-50.5,-51.0,-50.0,-46.0, 

*  -40.0,-34.0,-31.0,-32.0,-37.0,-37.5.-37.0.-30.0,-21.0, 

'  -21.0,-21,0) 

OATA  ( DL A T=  16. 0,14. 0,10. 0, 8. 0,6. 0.5. 0,4. 0,0. 0,-2. 0,0. 0,3. 0,6.0, 

'  la. a,  15. n, 16. oi 

OATA  (OLON= -110. 0,-1 06. 0,-1 05. 0,-10 4, 0,-10 2. 0.-98. 0,-98. 0,-103.0, 

►  -107.0,-109.0,-110.0,-113.0,-112.0,-111.0,-110.0) 


CALL  FOR  ENLARGED  PROJECTION  CENTERED  AT  WASHINGTON  ,  D.C.  AND 
A  PRIMARY  RANGE  ED'IAL  TO  20  00  N.  MILES.... 


CALL  3LOWUP  (  38  . 55, -77.  00, 2000. 0,2500.0, MOOTS, 

1  29HCENTERE0  AT  WASHINGTON  ,  D.C..29, 

►  30HPRIHARY  RANGE  EQUALS  2000  N.M.,30, 

►  18HPLOT  FOR  EXAMPLE  5,18, 

'  1,6. 0,6.0) 


C  COMPUTE  THE  SURFACE  AREAS  WITHIN  THE  2000  N.M.  RANGE  CIRCLE. ... 


CALL  SURFAREA  ( AL AT , AL DM , 19 , 3 0 . 0 , - 7 7 , 0 , NOP T S) 
CALL  SURFAREA  (OLAT.BLGN. 14,12,0, -95,0, MOOTS) 
CALL  SURFAREA  ( CL  AT , CLON, 29. 4 3 . 0 , -42 . 0 , NO FTS) 
CALL  SURFAREA  ( DL AT  ,  OLON , 1 5 . 6 , 0 , - 1 1 8 . 0 , NOP T 3 ) 

STOP 

END 


1 
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19 
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PRIMARY  RANGE  EQUflLS  2000  ND M 


PLOT  FOR  EXRMPLE  5 


CENTERED  AT  WASHINGTON  >  Do  Co 


FIG.  5:  FOUR  SURFACE  AREA  CALCULATIONS  WITHIN 
2,000  N.MI.  OF  WASHINGTON,  D.C. 
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- 

—  AREA  A  - 

SURFACE  AREA 

OF 

A 

EQUALS  290, 288,0 

SO. 

NAUTICAL 

PERIMETER 

OF 

A 

EQUALS  4,171.7 

NAUTICAL  MILES 

- 

—  APEA  8  - 

SURFACE  AREA 

OF 

3 

EQUALS  >16  1,348.4 

*  *  *  *  * 

SQ. 

NAUTICAL 

PERIMETER 

OF 

-3 

EQUALS  4,617.3 

NAUTICAL  MILES 

- 

--  AREA  C  - 

surface  area  of 

C< 

1) 

FOUALS  70,946.4 

*♦ *  ♦  *  ** 

SQ. 

NAUTICAL 

PERIMETER  OF 

C( 

1) 

EQUALS  1,151.7 

NAUTICAL  MILES 

SURFACE  AREA  OF 

C( 

2  ) 

EQUALS  640,o52.6 

*  *  *  *  *  *  * 

SQ. 

NAUTICAL 

PERIMETER  OF 

C( 

2) 

FOUALS  3,312.9 

NAUTICAL  MILES 

SURFACE  AREA  OF 

C( 

3) 

FOUALS  35,571.3 

♦  »  ♦  *  *  *  ♦ 

SQ. 

NAUTICAL 

PERIMETER  OF 

C< 

3) 

EQUALS  769.7  NAUTICAL 

MILES 

MILES 


MILES 


MILES 


MILES 


MILES 
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EXAMPLE  6 


The  triple  range  intersection  of  example  4  is  computed  here.  Since  the  area  of 
interest  is  small,  the  field-of-view  was  enlarged  for  better  observation  of  the  intersection. 
A  primary  range  was  not  required  and  was  set  equal  to  0.0. 

Triple  coverage  equals  17,995  square  n.mi. 
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PROGRAM  EXAMPLES  1 

DIMENSION  N0PTSC25) ,XLAT(15I ,XLON«15)  2 

3 

V 

C  STORE  COORDINATES  OF  INTERSECTION....  5 

6 

7 

OATA  (X  LA T=32.99, 12. 36. 31. 57, 30. 71, 29. ft  0,29. 83, 23. 69,28.29,29. 08,  6 

*  28. 32.20. 43.30. 59.il. 79. 32.99)  9 

DATA  (XL0N= 121. 95.122. 68,123.55,124. 32, 124.99,125.56,125.63,  10 

*  124.89,124.47,124.30,123.54,122.87,122.30,121.85)  11 

12 

13 

C  CALL  FOR  ENLARGED  PLOT  WITH  9L0M-UP  RANGE  EQUAL  TO  500  N.M .  14 

15 

16 

CALL  BLOWUP  (  29 . 0  0 , 1  23 . 50 , 0 . 0 , 5  0  0 . 0  ,  NOP  TS,  17 

*  27HSH0WING  TRIPLE  I HTERSEC T I  ON , 27 ,  18 

*  1H  ,1,  19 

*  18HPL0T  FOR  example  6,18,  20 

*  1,6.0,6.01  21 

22 

23 

C  COMPUTE  SURFACE  AREA  OF  INTERSECTION....  24 

25 

26 

CALL  SURFAREA  < XL AT , XLON , 1 4, 29. 00 , 1 24. 2 , NOPTS)  27 

29 

29 

STOP  30 

END  31 
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PLOT  FOR  EXRMPLE  6 


SHOWING  TRIPLE  INTERSECTION 

FIG.  6:  ENLARGED  TRIPLE  INTERSECTION  SHOWN  IN  EXAMPLE  (4) 
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******  ******* 


-  AREA  A 


9URCA0E  AREA 

OF 

A  Ff'UALS 

*  *  * 

17,994.4  99. 

*  *  *  * 

r<  Al.IT  I  CAL 

MILE" 

PERI METER 

Or 

A  F QUALE 

722.6  MAiJTICAL 

MT  L  c  9 

********  ***** 


-40- 


U 16275 


